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Abstract —We consider the uplink of non-cooperative multi¬ 
cellular systems deploying multiple antenna elements at the 
base stations (BS), covering both the cases of conventional 
and very large number of antennas. Given the inevitable pilot 
contamination and an arbitrary path-loss for each link, we 
address the impact of time variation of the channel due to the 
relative movement between users and BS antennas, which limits 
system’s performance even if the number antennas is increased, 
as shown. In particular, we propose an optimal linear receiver 
(OLR) maximizing the received signal-to-interference-plus-noise 
(SINR). Closed-form lower and upper bounds are derived as 
well as the deterministic equivalent of the OLR is obtained. 
Numerical results reveal the outperformance of the proposed 
OLR against known linear receivers, mostly in environments 
with high interference and certain user mobility, as well as 
that massive MIMO is preferable even in time-varying channel 
conditions. 

I. Introduction 

Multiple-input multiple-output (MIMO) systems have ap¬ 
peared as a method to increase the data rate and the spec¬ 
tral efficiency over wireless links through multiplexing gain 
and spatial diversity |[T]|. However, the evolution of cellular 
networks has resulted to a greedy need for higher data rates, 
achieved by means of multi-user MIMO (MU-MIMO) capable 
of exploiting multi-user diversity j2j. Despite that the channel 
state information at the transmitter (CSIT) contributes little to 
the multiplexing gain in single-user scenarios, its knowledge 
in MU-MIMO is crucial |3). 

Nevertheless, a new breakthrough technique, called massive 
MIMO or very large MIMO, makes it possible to apply low- 
complexity linear decoders to maximize system capacity by 
vanishing thermal noise, intra-cell interference, and estimation 
errors 0. In massive MIMO cellular networks, a base station 
(BS) deploys a large number of antennas, and multiple users 
can be served by the same time-frequency resource, since more 
degrees of freedom are provided. 

Unfortunately, the reuse of pilot sequences in adjacent cells 
causes pilot contamination §, which degrades the quality 
of CSIT. In addition, the commonly present mobility of the 
users, implying time-varying channel conditions, induces extra 
inaccuracy at the CSIT. As a consequence, the current CSIT 
cannot be perfectly known but only estimated. 

It has been shown that maximum-ratio combining (MRC), 
zero-forcing (ZF), and minimum mean-square error (MMSE) 
perform the same in single-cell systems, when the throughput 
per user is 1 bit per channel use, but MMSE is in general opti¬ 
mal |6j. The behavior of MMSE has been studied extensively 
in multi-cell networks, and in fact, in 0 an optimal linear re¬ 
ceiver (OLR), achieving the maximum signal-to-interference- 
plus-noise ratio (SINR), has been obtained. Nevertheless, some 
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limited prior work on massive MIMO has considered the 
effect of channel aging due to user mobility. Specifically, 
in [j8|, Truong et. al derived deterministic equivalents for 
the MRC (uplink) and maximum-ratio-transmission (MRT) 
for the downlink case, while in 0, 1 10J, the more complex 
deterministic equivalents for MMSE (uplink) and regularized 
ZF (downlink) were obtained and a comparison among the 
various strategies was provided. In the context of channel 
aging, the exact and approximate sum-rates have been obtained 
in |9]] for finite and large number of antennas, respectively, 
when zero-forcing (ZF) decoder is applied. The arising need 
is to derive an OLR that addresses the inevitable channel aging 
occurring in time-varying channels. 

Motivated by the above discussion, we hereafter derive the 
OLR in a cellular scenario subject to time-varying channel 
conditions, exploiting the correlation between the channel 
estimates and the interference from other cells. Furthermore, 
we elaborate on a comparison among various known receivers 
(MRC, ZF, and MMSE) to depict the outperformance of the 
proposed OLR and shed light on how the user mobility affects 
the performance in such case. In fact, we obtain lower and 
upper bounds of the sum-rate, accounting for channel aging, 
with an OLR assuming a finite number of BS antennas. In 
case of large number of antennas, we derive the deterministic 
equivalent approximation of the SINR that gives us the ability 
of avoiding lengthy Monte-Carlo simulations. Notably, it is 
of high importance to massive MIMO that in time-varying 
channels they are still favorable, but the performance of the 
system is limited, even if the number of antennas is cranked 
up, as shown by our results. 

II. System Model 

We consider a cellular uplink scenario, where the system is 
composed of L cells with one multi-antenna BS (A antennas) 
and K single-antenna users (A > K) per cell. All users share 
the same time-frequency resource, while we assume non- 
cooperative BSs. Taking into account a quasi-static block fad¬ 
ing model, where the channels are constant during the symbol 
period but vary from symbol to symbol, the channel vector 
< 7 iifc [n] £ C Nxl between the fcth user in the zth cell and the 
BS in cell l at the nth time-slot is decomposed as independent 
small-scale fading by means of huk £ C Nxl ~ CA/”(0 ,Ijv) 
and large-scale fading (shadow fading and path loss) by means 
of 3iik, given that the antennas are sufficiently well separated 
in realistic systems. More concretely, we have 

9uk[ n \ = y/lhikh Hk [n\. (1) 

In simple mathematical terms, the Ax 1 frequency-flat 
received signal vector at the (th BS (l = 1,2,... L ) at the ?rth 


time-slot is as follows 


L 

ViVA = y/P^Gii[n]Xi[n] + zi[n] 

L 

= Vp Hli l n ] D u 2x i M! (2) 

2=1 

where y/pXi[n\ £ C Kxl is the zero-mean stochastic transmit 
signal vector of K users allocated in the ?'th cell with p > 0 
being the average transmitted power per user, n; is the additive 
white Gaussian noise (AWGN) vector at the receiver (BS), 
modeled as CM (0, Ijv), and G u [n\ = [g H1 [n],... ,g HK [n]] £ 
C NxK denotes the complex combined matrix between all 
users in cell i and BS l. Note that in ((2}», we have applied <|T]) 
in matrix form, i.e., where Hu [n] £ C NxK is the fast 
fading matrix and Du is a K x K diagonal matrix with 
[Du\kk = Pli^\ 

The detection of the signals transmitted from I\ users by 
their associated BS necessitates the knowledge of channel state 
information (CSI), which can be obtained by using uplink 
pilots. Unfortunately, in realistic systems, the CSI at the BS 
is degraded by pilot contamination. This effect arises due to 
frequency reuse across a multiplicity of cells. According to 
MMSE channel estimation, the channel can be written as © 

Gu[n] = Gu[n\ + Gu[n\, (3) 

where vec^G^n]^ = Hu[n]Du ~ CM ^0, Ijy ® £>u^ and 

vec (Gljn])~CAA(0 , Ijv 8 (Du — -^»)) are the indepen¬ 
dent channel estimate and channel estimation error matrices, 
respectively. Note that Hu [n] is a N x N matrix with Gaussian 
entries having zero mean and unit variance as well as Du is 
a K x K diagonal matrix with elements [Du\kk = Puk = 

Q 2 

—, where p p is the transmit power of each pilot 

symbol (p p = rp). Nevertheless, an important relationship, 
going to help our derivations, concerns the estimated channel 
between the /th BS and the users in cell i. It is given by 0 

Gu[n] = Gu[n]Ru, (4) 


where G h [n} = \g H1 [n ],... ,g liK [n]] £ C NxK is the esti¬ 
mated combined channel matrix from all users in cell i to BS 
l and K„ =diag{|ff,|g,..., |fj}. 

Besides pilot contamination, in any common propagation 
scenario, a relative movement takes place between the anten¬ 
nas and the scatterers that degrades more channel’s perfor¬ 
mance. Under these circumstances, the channel is time-varying 
and can be modeled by the famous Gauss-Markov block 
fading model 03’ which is basically an autoregressive model 
of certain order that incorporates two-dimensional isotropic 
scattering (Jakes model). More specifically, our analysis targets 
to relate the current channel state with its past samples. 
For reasons of simplicity and computational complexity, we 
consider the following autoregressive model of order 1 J8) 


Gu[n\ = olGu[ n - 1] + E u [n], 


(5) 


where vec^EjJn]^ ~CM (0, Ty ® (l —a 2 ) Du) and Gu[n— 1] 
are uncorrelated. They model the stationary Gaussian channel 
error vector due to the time variation of the channel and the 
channel at the previous symbol duration. The a = Jo (27 t foT s ) 
parameter, where Jo(-) is the zeroth-order Bessel function of 
the first kind, fjj and T s are the maximum Doppler shift and 
the channel sampling period, expresses the two-dimensional 
isotropic scattering. It is worthwhile to mention that the 
maximum Doppler shift fp equals Jd = , where v (in 

m/s) is the relative velocity of the user, c = 3 x 10 8 m/s is the 
speed of light, and f c is the carrier frequency. 

The combination of 0 in 0 allows the characterization 
of both effects at the same time instance, i.e., it provides the 
channel at time-slot n, accounting for pilot contamination and 
time variation of the channel. Specifically, we have 

Gu [n] = uGu [n - 1] + E u [n] 

= aGu [n - 1] + Eu [n], (6) 


where Gu[n - 1] and Eu\n] = aGu[n - 1] + Eu[n] ~ 
CM ^0, Ijv <8 ^Du — a 2 Duj) are mutually independent. 

We focus on the study of linear detection by the receiver 
(BS) by means of a linear N x K matrix Wj[n] dependent on 
the channel estimate Gu[n] After applying the detector W;[n] 
to its received signal y I [n] at time-slot n, the /th BS obtains 

L 

n = W'l [n\y l =^W“ [n] G u [n\x, [n] + W" [n]zi [n] 

i—1 

L 

= a y /pW H l [n\'^ j Gu[n - l]«»[n] 

2=1 

L 

+ VpW'i' [n] y; Eli[n]xj [n] + W'lWzPn] (7) 

2=1 

where in ([7]), we have used 0- Hence, the kth element of r/, 
corresponding to the signal from the kth user, is given by 

rik = aVpw H i k [n]gn k [n - 1 \xi k [n] 

K L 

+aVpY}"ik[n\9iijl n ~ l ] x ij W\+ a VpYl w ^^ li1 ] Xi W 

jy£k ly£l 

L 

+ VpY 1 w "k WEu [n]xi [n\ + w" lk Zi [n], 

2=1 


where xi k and Wn- are the A th element and column of Xi 
and Wu respectively. Taking into account the independence 
between Gu[n — 1] and Eu[n], the achievable uplink SINR 
SINRfc of the A'th user can be written as in ([8][_] Note that 

Rp, = i.e., 

i-l 


L L K 

Re = dV- n = 53 \fi lik ~ offilik 


l N - 


2=1 k—1 


'Regarding shadowing, i.e., ftifc, it is independent of symbol’s index n 
because we assume that the coherence time is large relative to any large-scale 
delay constraint of the channel. 


2 We assume encoding of the message over many realizations of all sources 
of randomness in the model including noise and channel estimate error jTJ. 




SINR fc = 


a 2 p\w H lk [n\g nk [n - l]\' 


a2 PJ2f^k KfcWfliZjt” - l )\ 2 + PEt=i w "ki n } R E w ik[n} + u 2 PY,Ui II wH iM G u[n - 1]|| 2 + I K fe M 


( 8 ) 


III. Optimal Linear Receiver 

In this section, following a similar procedure as in |7), we 
derive the OLR maximizing the received SINR. Specifically, 
we have 


SINR fc 


l^fc[n]g;; fc [n-l]| 2 

w H lk [n]^ k [n-l]wi k [n] 


w"i k [n]^ 2 [n-l} 

2 

s *2 /2 [n-i]g H „.[n-i] 

2 

™iki n 

‘—‘/c 

n—l\wi k [n] 


= 9ak[n ~ ^^[n-llgukin - 1], 


(9) 

( 10 ) 


where we have denoted that Efc[n — 1] = 

E & k 9iij[n - l]9uj[n ~ !] + EjyGji[n - ^[n - 
!] + 5? (E<=i hi + Xj I,y as well as in ([9]) we have used 
Cauchy-Schartz’s inequality. According to its property, the 
equality holds when Wik[n] = cEj^ 1 [n — l\g ak [n — 1] for any 
nonzero c £ C. It is important to mention that Efc[n — 1] can 
be also written as 


L 

Sfc[n— 1] =5^Gti[fc] [n-l}G" li[k] [n- 1] +<7 2 Ijv, (11) 

i =1 


Application of the eigenvalues decomposition to the follow¬ 
ing N x N matrix 


L 

S[n ~ !] = l n “ \ n ~ !] 

i =1 

— Gii[k][n - 1 ]Ri{k]Gii[k] l n ~ 1] 

- Him l n ~ \ n - *] 

with R l[k] = J2i=i R u{k] and D i[k] = Ei=i Aiffc] provides 
the eigenvalues in terms of a diagonal matrix B[n — 1] and 
the corresponding eigenvectors by means of the columns of the 
unitary matrix U as S[n — 1] = U H [n — 1 ]B[n — l]U[n — 1], 
Given that the rank of S\n — 1] is K — 1, the matrix B has 
the form 

N-K+l 

B[n- 1] = diag[Ai,..., X K -i, 0,..., 0 ] [n _ 1] , 

where the first K —1 terms represent the respective eigenvalues 
of S[n — 1] and the subscript [n — 1] in front of a matrix 
denotes its instance at the (n — l)th slot. Taking advantage of 
the eigenvalues decomposition and ( p~T| ), we can express © 
as a sum of two independent terms as follows 


where Gu\k] [n— 1] is the matrix Gu [n— 1] with its fcth column 


removed and a 2 = 


a 2 (Ej=lGi 


Z,= E 


E ll E ll 


+^E[G h G^] 

i^l 

K ( L \ 

= ^ ( ^ fink I In, 

) 

and the expression of the typical MRC is Wi k [n]=ag uk [n—l], 
the SINR for the OLR and MMSE are expected to behave 
similarly in low interference conditions, while the proposed 
OLR achieves higher sum-rate in interference-limited scenar¬ 
ios, as expected |7J. However, the main contribution rests on 
the exhibition of the effect of the user mobility on the OLR, 
which is clearly described in ©■ As can be seen, the SINRs 
corresponding to all receivers under investigation coincide 
after a certain value of a that makes the term multiplied by 
1/a 2 negligible. 


SINRfc = (u[n— 1 ]g Uk [n- 1]) " B[n- 1] + 


Remark 1: Given that the typical MMSE receiver is 


Wi k [n]=a[Gu[n-l}G H ll [n-l]+ “2 + Ijv) g llk [n-i\, 


where the matrix Zi is deterministic and equals to 

L 


x U[n - l]g uk [n - l] 

= [ftV > ■ • • > ft* k N ][„_!] 

9llki 9llkK—i 9llki, 


gnki\ 


1 T 


Ai + O' 2 


K-l i- 


Xk-i + ct 2 ’ <t 2 ’ ’ a 2 


[n-l] 


N-K+l 


_ y' \gukj\ 2 i | 


1=1 


lk+ 


( 12 ) 


i=i 


where U[n - 1 }g llk [n - 1] = [g ukl ,... ,g UkN j^,. The 
multiplication of gu k [n — 1] with the unitary matrix 17[n — 1] 
preserves the statistical properties of the former one. Moreover, 
A j for j £ [1, K — 1] and \gn k | are independent, since Xj 

represents the eigenvalues of S[n— 1], i.e., Gu[n — l ]G h [n- 1] 
with the fcth column removed. 

The sum-rate R k of the fcth user is written as 


R k = E [log 2 (1 + SINRfc)]. 


(13) 


Because of the arising difficulty to obtain the exact distribution 
of the SINRfc, we resort to derive bounds of the sum-rate that 
prove to be very tight, as shown. 



























A. Upper Bound 

Proposition 1: The upper bound of Ur for user k in the 
reference cell j with OLR, accounting for user mobility, is 

K-1K-1 
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yN-K+u-r 1 ^r+K-N-2 


Substitution of ( fl9| and ( [20] ) into ( fl5| concludes the proof. ■ 
B. Lower Bound 

Proposition 2: The lower bound of lit, for user k in the 
reference cell j with OLR, accounting for user mobility, is 
given by 


Rk > log 2 ( 1+2 (K — 1 ) exp 


-27- 


1 


211 “ *Ii) 


K -1 if-1 n .K-N-2 
< yg 1J lvuT' l% 


(. N-K+u)\ 


T(N -K + u + 1) 


(14) 

where Ei (z) = — f x e~ t /tdt is the exponential integral 
function | |T 6 ] Eq. (06.35.07.0002.01)]. In addition, ti v is the 
rith element of the diagonal matrix Dim and Di vu is the 
(u, u)th cofactor of a (K — 1) x (K — 1) matrix Di whose 
(i,j) th entry is {D+j = t 3 u l . 

Proof: From ( | 1 3j ), we have 

Rk < log 2 (1 + E [SINRfe]) 

, ll 0 g 2 (l+E[I 1 +I 2 ]), (15) 

where we have applied Jensen’s inequality, since log 2 (•) is a 
concave function. 

The expectation in GU can be calculated by taking the 
expectation of each term separately. Hence, the first term can 
be written as 

Epr,] = (K — 1 ) E 1 ^ 




2 \ N—K+u +1 


N-K+u 


gt;„ _|_ Ei r + l( — 

r=0 ' lv 


( 21 ) 


where 7 is Euler’s constant with numerical value 7 ~ 
0.577216 and Ei„ (z) = e~ zt /t n dt, n = 0,1,2,..., 
Re( 2 ) > 0, is the exponential integral function of order 
n p 6 j Eq. (06.34.02.0001.01)], which can be related to the 
incomplete gamma function (Ei„ (z) = x n_1 r (1 — n,x)j. 

Proof: The proof relies on the general bounding technique 
according to which we can re-express (p~3|), according to 


Rk = E 


> E 


log 2 (l + exp (in (Ji + J 2 ))) ] 

log 2 (i + exp^(In (210 + ln ( 2 ^ 2 ))) / 2 ^ 


> log 2 [l + exp(E[(In (210 + In ( 21 2 ))] / 2 ) \ ( 22 ) 


06) 

(17) 


where first we have used the arithmetic mean-geometric 
mean inequality, and then we have exploited the fact that 
log 2 (1 + exp (x)) is convex in x and thereafter applying 
Jensens inequality. The last expression is constituted by two 
parts. Regarding the first part, we have 


k- 1 


E 


(18) 


[ln(2J0>E[E ln 


1=1 


2 (K-l)\g llk . 

A j + a 2 


'{K- 1) 


(23) 


= E 
= E 


In 


2 (K - 1 ) | g x 


Ik | 


A + O 2 

In (2 (K- 1) \guk\ 2 ) -E[ln(A + < 7 2 ) 

r°° i 

/ ln ( 2 ( K ~ 1 ) a:e x ) dx - 
Jo 


n ,< 7 c) 


(19) 


K-1K-1 „ .K-N-2 

x U lvu*i v 


where in we have used that the terms of the sum are 
i.i.d. Given the independence between \guk\ and A and that 
\giik\ 2 is an exponential variable with unit mean and unit 
variance, we lead to ( fl7] >. Keep in mind that the K— 1 nonzero 
eigenvalues of S[n -"Tj have the same distribution as the N 
eigenvalues of S H [n — 1], Specifically, substituting the PDF 
of the unordered eigenvalue of S H [n — 1 ], obtained by [ 13 1 
into (fl7]l, ( p~ 8 ] > is obtained. The integral in ( fT 8 ] is solved by 
using |(15] Eq. (3.353.5)]. Concerning the second term in 
we have, following the same way, that 

E [J 2 ] = (iV ~f ~ 1} E [\g llk \ 2 ] = (20) 


^ ^ T(N-K + u + l) 

f ln (A + a 2 ) X N ~ K+U e~ nk dA 
Jo 


= ln (2 (K— 1)) — 2 7 
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n 


K -1 
i<j 


K- 1 K -1 

EE 


n . K-N-2 

rJlvu b i v 


(24) 

(N-K+u)\ 


(*«-*«) m-K+u+i) 

2\ 9 N—K+u / o 

O' \ i\t _ +— X ' ( O 


Note that in ([23]>, we have used again the arithmetic mean- 


geometric mean inequality and in (|24|) the fact that \gi 


Ik 



































is exponential with unit mean as well as that A is an 
unordered eigenvalue of semi-correlated Wishart distribution 
according to m- The last expression is obtained after us¬ 
ing ['15, Eq. (4.332.1), (3.351.3)] and (T4) Eq. (47], Similarly, 
E[ln(X 2 )] is obtained as 

E[ln (2T 2 ) ] = -2 7 + In (2 (K - 1)), 

which when substituted with ( [25] ) into ( | 22 ] >, returns ( | 2 Tj ). ■ 

IV. Deterministic Equivalent of SINR 

The emerging technology of very large MIMO leads us to 
study the effect of finite Doppler shift on their performance. 
Specifically, the following proposition quantifies the degrada¬ 
tion of the SINR due to the relative movement of users. 

Proposition 3: The deterministic equivalent of SINR/,, for 
user k in the reference cell j with OLR, accounting for user 
mobility, is given by 

SINRfcX ti'TuPuk, (26) 


where 


-l 




In, 


and 5nk = lim^oo 6)^, where for t = 1 , 2 ,... 
K A i / L 


C = E 


\j =i 



PukN 


with initial values 5^ = a 2 / (j2i=i &U + ), 
Proof: It is obtained by means of ([TO} 


AJ for all k. 


SINR fc = g" llk [n - l]S fc 1 [n]g llk [n - 1] 


1 [tt,— T]/3 zzfe 


:trT 


liPllk, 


(27) 

(28) 


where the symbol x relates two infinite sequences and is 

equivalent to a,S ’ - > denoting almost sure convergence. As 

N—fOO _ _ 

far as the derivation is concerned, ( ]27| i follows from fl7| 
Lemma B.26] as well as in ( |28| > we have applied the rank- 
1 perturbation lemma Jl 8 | | and then | [I9[ Theorem 1], ■ 

The deterministic equivalent uplink sum-rate can be ob¬ 
tained by means of the dominated convergence [ 20 ] and the 
continuous mapping theorem pT| as 

Rk - log 2 (1 + SlNR fc ) ———A 0. 

N—foo 

V. Numerical Results 


This section provides the representation of the impact of the 
user mobility on the sum-rate as well as the assessment of the 
proposed bounds and deterministic equivalent. The scenario 
under consideration includes L = 7 cells with K = 10 users 
per cell. For reasons of comparison, the large-scale fading 
matrix Du,,i = 1,7 and the rest parameters, are chosen 


|7£ Specifically, T = 196, r = 10, p p = p, and 


SNR = p. Taking as reference cell 1, the behavior of the 


3 This large-scale matrix models practical channel conditions in hexag- 
onal cells with randomly distributed users under log-normal shadowing 
(lnA/*(0, 8)) and path loss with exponent equal to 4. 


spectral efficiency in bits/ sjHz is investigated. In particular, 
we study the total spectral efficiency given by 

K 

k =1 

where R k is the sum-rate of user k in cell 1 and T is the 
coherence time interval (in symbols). Note that Du is given 
in terms of (3 expressing the effect of the interference from 
other cells. It has to be stressed that in all cases the outperfor- 
mance of the OLR is depicted. In addition, the deterministic 
approximation of the sum-rate is indistinguishable comparing 
to the simulation results. It is worthwhile to mention that the 
sum-rate saturates for high SNR values (interference-limited 
system), i.e., the system performance cannot be improved by 
using more power, since increase of the transmit power results 
also to increase of the interference from other cells. 

In Fig. 1, the simulated achievable sum-rate of the OLR 
along with its deterministic equivalent ( |26| ) are plotted against 
the average SNR p. In addition, the typical MRC and MMSE 
decoders are simulated as well. Addressing different inter¬ 
ference configurations across different cells, i.e., [3 = 1 and 
[3 = 4, we observe that the OLR presents the best performance, 
which increases as shown by the zoomed gap, when the 
intercell interference increases. The explanation rests on the 
fact that the correlation between intercell interference and the 
channel estimate are considered by the structure of the OLR. 

In Fig. 2, the effect of the user mobility by means of 
the Doppler shift fdT s on various receivers is illustrated by 
depicting the variation of the sum-rate with the normalized 
Doppler shift. It can be observed that the Doppler shift 
extensively limits the performance of the channel. In particular, 
the sum-rate decreases from almost 9.2 bits/s/Hz to zero 
by following the form of the Jo (•), while the Doppler shift 
increases. Moreover, we observe that the gap between the 
performances of the various receivers vanishes to zero for 
values of a greater than 0.25. Nevertheless, increasing the 
number of antennas and keeping the other parameters fixed, the 
sum-rate becomes higher for all values of JdT s , as expected. 

The tightness of the proposed bounds is demonstrated in 
Fig. 3, where iQ 4), gg, and ( |26l > are compared with the 
simulation results. We can easily observe that the proposed 
approximations are very tight across varying SNR. 

VI. Conclusion 

We have proposed an OLR capable of accounting for 
not only the detrimental effects of pilot contamination and 
path-loss, but also the critical time variation of the channel 
due to the movement of the users. Taking into account the 
correlation between the channel estimate and the interference, 
we showed the outperformance of the OLR with comparison 
to MRC, ZF, and MMSE by deriving closed-form lower and 
upper bounds of the sum-rate for finite N. Moreover, we 
obtained a deterministic equivalent of the SINR when both 
the numbers of users and antennas go to infinity with a given 
ratio under the same channel conditions. Simulations witness 
the tightness of the bounds as well as the validity of the results. 
It is worthwhile to mention that the deterministic equivalent 
provides a tight approximation of the SINR, allowing us to 
avoid unbearable time consuming Monte Carlo simulations 
when the number of BS antennas is very large. 











Fig. 1. Sum spectral efficiency versus SNR for different interference 
configurations (/3 = 1 and (3 = 4). 
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Fig. 2. Sum spectral efficiency versus normalized Doppler shift for different 
numbers of BS antennas ( N = 50 and N = 100). 
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